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SAISS  XW  NUCUBAfi  SEACTOB  COMPONSirrS 


'  Llojrd  G.  Alexander 
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AB3T1AGT 

Knowledge  of  local  heating  rates  is  needed  for  estimation  of 
operating  tes^ratures^  thermal  stresses,  and  cooling  requirements 
in  nuclear  reactor  coi^)onents.  Heat  is  liberated  bgr  the  dissipation 
of  the  energjr  of  fission  fragments,  beta  particles,  fast  neutrons, 
and  gamma  photons*  Heating  rates  are  formulated  in  terms  of  either 
neutron  or  gamma  flux  densities,  the  corresponding  coUisicui  probabilities, 
and  appropriate  ^energy  transfer  functions,  the  forms  of  idiich  are  deriyed. 
Special  methods  of  estimating  the  flux  densities  are  discussed.  The  data 
on  the  magnitudes  of  the  yarious  energy  sources  are  reyieved* 


OM  THE  CALCULATION  OF  THE  DISTRIBUTION  OF  HEAT  LIBEHATION 
RATES  IH  MUCIEAR  REACTOR  COMPONEHT£ 


By 

Lloyd  G»  Alexander 
Oak  Ridge  National  Laboratory 


INTRODUCTION 

The  designer  of  nuclear  reactors  needs  to  knov  the  local  rates  of 
heat  liberation  in  the  various  reactor  components  in  order  to  estimate 
temperatures,  cooling  requirements,  thermal  stresses,  corrosion  rates, 
radiation  damage  annealing  rates,  tenqperature  dependant  reactivity- 
changes,  a  nd  so  ono  The  heat  liberation  rates  can  be  formulated  in 
terms  of  the  local  neutron  and  gamma  fluxes^  ^jj-and  Admittedly, 

values  of  these  are  not  easily  come  by,  Ne-vertheless,  recent  progress 
in  neutron  and  gamma  particle  transport  theories  makes  it  possible  to 
treat  many  systems  of  engineering  interest.  It  is  the  purpose  of  this 
paper  to  formulate  the  beat  liberation  rates  and  to  indicate  methods 
of  estimating  the  fluxes « 

EtffiRGY  SOURCES  AND  MODES  OF  LOCAL  DEPOSITION  OF  HEAT 

Heat  is  liberated  in  reactor  components  by  the  dissipation  of  the 
kinetic  energy  of  fission  fragments,  beta  particles,  gamma  photons,  and 
neutrons.  These  particles,  in  penetrating  matter,  collide  with  elec¬ 
trons  «nri  atomic  nuclei  and  impart  energy  to  them.  In  same  collosion 
processes,  part  of  this  energy  is  reradiated,  but  for  the  most  part  it 
is  rather  rapidly  degraded  into  thermal  motion  of  nuclei  and  electrons 
in  the  macros  copies  l].y-near  vicinity  of  the  collision. 
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Fission  fragments,  primary  gamma  photons,  and  neutrons  result 
directly  frcm  fission.  Fission  fragments  emit  beta  particles  and 
gammn  photons  during  radioactive  decay.  Heutrons  diffuse  through 
the  reactor  and  generate  gamma  photons  by  inelastic  c  ollisions  with, 
and  absorptions  in  various  nuclei.  Such  activated  nuclei  often  decay, 
emitting  beta  and  gamma  particles.  Photons  arising  from  all  these 
processes  may  experience  many  scattering  collisions,  losing  part 
of  their  energy  at  each  collision, 

FISSION  FRAGMENT  AND  BETA  DECAY  ENERGY 

According  to  Shapiro  (  16),  tw  and  occasionally  three  heavy 
particles  result  from  the  fission  of  by  thermal  neutrons,  Hanna 
(  9  )  critically  examined  the  data  as  of  1950  and  concluded  that  the 

most  probable  value  for  the  average  kinetic  energy  of  the  fission 
fragments,  lies  in  the  range  168  ^  5  mev  per  fission. 

Because  of  their  high  charge,  fission  fragments  have  short  ranges. 
(  16  ).  ..Accordingly,  their  energy^  Gffjls  dissipated,  macroscopically 
at  the  point  of  fission, 

K.  Way  (  21  )  estimated  the  e  i:|*gy  released  by  the  beta  decay  of 
fission  fragments  ,%f^5  to  be  7  mev  per  fission.  Effective  ranges 
are  of  the  order  of  1  mm,  or  less  in  reactor  materials j  hence,  the 
beta  energy  is  also  liberated  very  near  the  point  of  fission. 

The  sum  of  the  fission  fragment  and  beta  ray  energies  is  thus 
about  175  mev  per  fission.  The  heat  release  rate  is  proportional 
to  the  nvimber  of  fissions  per  unit  volume  per  unit  time  and  may  be 
formulated  as  follows; 
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X  d  E  mev. 

Here,  (E«) represents  the  energy  dependent  neutron  flux —  l.e., 
total  path  length  (cm)  traversed  ty  neutrons  having  energy  E,^(iaeT)  per 
unit  volume  (cm.^)  per  unit  time  (sec,)  per  unit  energy  range  (mevj  , 

The  integration  is  carried  over  all  energies;  however,  there  are  negligibly 
few  neutrons  having  energies  less  thaji  that  corresponding  to  the  tempera¬ 
ture  of  the  medium  (30°  C.  -v-- 0.025  ev,),  or  greater  than  10  mev. 

In  •^thermal"  reactors,  the  great  majority  of  the  flux  is  concentrated 
into  a  narrow  band  around  the  temperat\ire  of  the  medium.  Furthermore, 
the  fission  cross  sect  ion,  decreases  rapidly  with  increasing  neu¬ 

tron  energy.  Accordingly,  in  thermal  reactors,  the  integral  is  appro¬ 
ximated  very  well  ty  the  product  of  the  “thermal"  flux  (say  the  flux 
lying  between  O.ol^^d  0.1  mev)  and  the  fission  cross  section  evaluated 
at  the  temperature  of  the  medim.  I^t  the  thermal  flux  be  denoted  by 
(th),  whence  /^0.1 

(th)  -  /  (In)  dEn  ^ 

‘^0.025 

The  product  of  the  first  five  qu^ities  of  Equation  1  is  commonly  denoted 
W  Zf  (En)j  the ‘Wcroscopic  fission  cross  section".  It  is  the 

probability,  per  cm.  of  travel,  that  a  neutron  of  energy  E  will  experience 
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a  fission  collision,  and  has  units  of  fissions  per  neutron-cm.  Using 
Equation  2,  and  converting  to  BTU/hr.ft.^,  Equation  1  takes  the 
form,  for  thermal  reactors 

Qf  -  (Off  Sff^)  4  (ti>)  3 

For  U^^5  equals  549.  Taking  m  to  be  1  gm.  U235/  cc,  gives  ot 
r,  of  1,40  fissions/neutron-cm.  At  a  thermal  flux  level,  J?^(th),  of 
10^3  neutron-cri./cm.3  sec,,  the.  heat  production  is  about  400  watts/cc 
or  40  X  10^  BrU/hr,ft.^ 

REUTRON  ENERGY 

Neutrons  Impart  kinetic  energy  to  atomic  nuclei  by  three  processes; 
(a)  absorption,  (b)  elastic  scattering,  (c)  inelastic  scattering. 

In  an  absorption  collision,  the  neutron  enters  and  combines  vd.th  tl» 
target  nucleus,  forming  the  next  higher  isotope.  The  new  nucleus  is 
"knocked  on",  and  its  velocity  is  readily  calculated  from  a  simple  momen¬ 
tum  balance. 


The  thermal  motion  of  the  nucleus  prior  to  the  collision  may  be  neglected, 
since  it  is  very  small  compared  to  the  neutron  velocities  of  interest. 


Therefore, 

> 

mVns 

(M+m;  Vjj/ 
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where 

m  - 

atomic  mass  of  neutron  (m«rl) 

M  - 

atomic  mass  of  nucleus 
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Vn  =  velocity  of  neutron 

velocity  of  nucleus  after  collision 
^y  definition,  the  kinetic  energies  of  the  nucleus  and  neutron  are, 
respectively 


%i  e  i  (m-in)  Vji? 

5a 

5b 

Combining  Equations  4»and>5 

—  mBn  /(M+m) 

6 

Defining  S^c  fraction  of  the  kinetic 

energy  of  the  incident 

neutron  appearing  as  kinetic  energy  of  the  target  nucleus  after  a  cap¬ 
ture  collision,  one  has 


Snc  -  =  m/(Mf  m)  7 

The  remainder  of  the  neutron  kinetic  energy  plus  the  binding  energy 
is  released  in  the  form  of  a  series  of  so-called  “capture*  gamma  photons 
of  decreasing  energy.  The  N'  isotope  may  be  radioactive,  and  may  decay 
by  alpha,  beta,  and  gamma  emission.  Later  sections  deal  with  heat  re¬ 
lease  from  these  sources. 

The  local  heating  effect  due  to  the  capture  of  fast  neutrons  is 

thus 


vOD 


^  /  Snc  ^nc  « 


'Bo 


5 


where  is  the  so-called  macroscopic  neutron-capture  cross  section.  It 
is  the  probability,  per  unit  length  of  travel,  that  a  neutron  will  ex¬ 
perience  a  capture  collision.  The  lower  limit  of  integration,  Eq,  is 
chosen  suffiently  high  that  Equation  8,  derived  on  the  assumption  that 
the  energy  of  the  target  nucleus  prior  to  the  collision  is  negligible 
in  comparison  with  the  initial  neutron  energy  En,  remains  valid,  but 
also  sufficiently  low  that  only  a  negligible  portion  of  the  neutron  energy 
remains  to  be  transferred  at  energies  less  than  Eq.  The  choice  depends 
somewhat  on  the  spectral  distribution  of  the  fast  neutron  flux,  but 
values  between  10  and  100  ev  are  nearly  always  satisfactory. 

If  two  or  more  kinds  of  atoms  are  present  in  the  scattering  medium, 
the  corresponding  Gnc  must  be  summed  over  all  of  these.  This  is  possible 
because  the  capture  cross  sections,  &re  additive. 

Both  momentum  and  kinetic  energy  are  conserved  in  elastic  collisions. 
The  appropriate  balances  are  most  conveniently  formulated  in  the  center-of- 
mass  system  of  coordinates.  By  combining  these  and  transforming  the  result 
to  the  laboratory  system  of  coordinates,  Glasstone  and  Edlund  (5)  obtain  a 
result  that  may  be  put  into  the  form 

E  ’  /E  =  m2  +  2tfoiGo3  6  m^  9 

^n/^n  ■(m'+  mj2 

where  0  is  the  angle  of  scattering  in  the  center -of-mass  system.  Since 
kinetic  energy  is  conserved  in  the  collision,  it  follows  that 

En/^^  =  1  -  10 

The  scattering  is  isotropic  in  the  center-of-mass  system,  hence  the 
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probability  P(^)d-^  that  a  neutron  will  be  scattered  elastically 
through  an  angle  ^  Into  angle  d  Is  related  (  5  )  to  the  angle 


^  hy 


p(-^)d-^  -  Sln^d^/2 
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By  definition,  the  average  fraction  one  of  neutron  energy  En  con- 
verted  into  target  nucleus  kinetic  energy^  is  given  t?y  the  integral 

Sne  =  /  (1  -;^n)  H^)  ^2 

l/o 

which  leads  to  the  simple'  res\ilt  that 
^ne  =  (l-oC)/2 

where  oC  -  -  m)  /(M  i-m)}  2 

The  elastic  collision  heating  rate  is  fomrulated  similarly  to  Equation 
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13. 

One  = 

The  lower  limit  Eq  again  lies  in  the  range  10  to  100  ev.  As  with 
Gjjjg  must  be  summed  with  respect  to  all  kinds  of  atoms  present. 

If  the  compound  nucleus  formed  hy  the  interaction  of  a  neutron 
with  a  target  nucleus  decomposes  by  the  emission  of  another  neutron, 
the  net  result  is  equivalent  to  the  scattering  of  the  incident  neutron. 

The  nucleus,  however,  may  bo  left  in  an  excited  state,  in  which  case  kinetic 
energy  is  not  conserved  and  the  scattering  is  Inelastic,  There  is  thep^ 
for  this  case,  no  relation  between  the  scattering  angle  and  the  energies 
of  the  incident  and  scattered  neutrons.  Instead,  with  heavy  (M  ^40  ), 
non-magic  nuclei,  the  probability  P(En,E^a)dE^  that  a  neutron  of 
energy  Ejj  will  be  scattered  through  an  energy  range  (E^  -Ej^)  into 
energy  range  dE’^  by  a  nucleus  of  atomic  number  Z^is  given  by  the 
Welskopf  (  7  )  formula,  which  may  be  put  into  the  form 


J  ^e  Eji  ^ne  (En)  Pa  (En)dEn 

'-'E- 
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P(EniEA,Z)dEi  =  fep/T^)  e3^(-EyT^) 

1  -  (1+En/T2;)  exp(-En/T2) 

where  is  the  so-called  "nuclear  temperature"  which  must  be  determined 
experimentally  for  each  nucleus.  Values  are  listed  by  Bjorklund  (2)  idio 
reports  that  T2  is  independent  of  and  and  is  a  function  of  only  Z, 
the  atomic  number. 

The  scattering  is  isotropic  in  the  center -of -mass  system  of  coordinates 
However,  with  heavy  nuclei,  it  is  also  approximately  isotropic  in  the  labora- 
toiy  system.  Accordingly,  the  probability  P(9)  dO  that  the  neutron  will  be 
scatted  through  angle  9  into  angle  d9  is  given  by  Equation  11,  where  9  is 
now  the  scattering  angle  in  the  laboratory  system.  Specifying  E^,  E^k  and 
9  completely  determines  Er,  the  kinetic  energy  of  the  target  nucleus  (for 
the  case  of  negligible  initial  energy) .  A  simple  momentum  balance  readily 
yields  the  relation 

*  * 

Ek/Bq  =  (1/M)  [1-2  (E^)^/^^  Cos  0  +  (EyEn^  17 

By  definition,  the  fraction  <5^ni  incident  neutron  energy  E^  that 
appears  as  kinetic  energy  Ejj  of  the  target  nucleus  is  given  by 

P - - - ^ 

,(%/E^)P(E„,E-^I^(En,Z)  dEi  P(0)de  18 

? 

where  2ni(^»^)  inelastic  scattering  "cross  section". 

Substitution  from  Equations  11,  I6,  and  17  followed  by  integration  gives 

r^i(En,Z)  =  2:ni(En»Z)Tz  -  (3+2Tz/E^+2En/rz)  exp(-E^/T22l  19 

M  Ql  -  (1+En/Tz)  exp  (-Eu/T2)J 


,S'ni(En-Z)  = 
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Equation  19  is  valid  only  for  heavy,  non-magic  nucleic  For 
light  and  magic  nuclei,  the  possible  energy  levels  in  the  excited 
target  nucleus  aire  more  or  less  discrete,  and  accordingly  the  pro¬ 
bability  P(Ejj,Ej^)  is  distributed  in  discrete  bainds  which  must  be  deter¬ 
mined  experimentally  for  each  nucleus  (7  )«  The  data  are  then 

applied  to  Equation  19  in  place  of  Equation  l6o 

For  all  nuclei,  5nl(®n>2)  is  a  function  of  Ej^o  The  local  heating 
rate  due  to  inelastic  scattering  of  fast  neutrons  accordingly  takes  the 
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where  Eq  is  the  threshold  ener^  for  inelastic  scattering.  must  be 
summed  with  respect  to  all  kinds  of  atomic  nuclei  present. 

ACTIVATED  NUCLEI  BETA  DECAY  EMERGY 


Nuclei  rendered  unstable  by  the  absorption  of  a  neutron  may  emit 
alpha,  beta,  or  other  particles  as  well  as  gamma  photons.  Beta  emission 
is  the  most  common  mode  of  decay.  For  a  given  transition,  the  beta 
particles  are  emitted  with  energies  ranging  continuously  from  zero  up 
to  some  characteristic  maximum,  (max),  which  is  the  most  commonly 
tabulated  value  (  10  ),  Since  the  contribution  to  the  local  heating 
from  beta  decay  of  activated  nuclei  is  usually  small  compared  to  other 
contributions,  it  is  generally  sufficiently  precise  (  15  )  to  take  the 
average  electron  energy  to  be  one  third  (max).  Thus,  the  local 
heating  rate  due  to  ”  neutron-captur^' betas  becomes 


*^ncy7  " 


X^nc^^n^^n't^u)  dEn  ^nJpi(max)/3 
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7  where  ni  is  the  nia^r  of  betas  cf  characteristic  energy  Eei  (max)  emitted 

y  A+K 

in  the  decay  per  neutron  capture.  Since  isomeric  transitions  occur, 

A 

the  ni  are  numbers  equal  to  1  or  less.  The  decay  process  may  proceed 

-f  /te 

throu^  several  atomic  nianbers.  Note  that  ^contribution  from  ”thennal" 
neutrons  in  Equation  21  may  be  large  and  override  the  contribution  from 
fast  neutrons. 

GAMM4  PHOTON  ENERGY 

Gamma  photons  penetrating  matter  exhibit  three  interaction  processes; 
(a)  photoelectric  effect,  (b)  pair  production,  and  (c)  Con^jton  scattering. 
In  the  first  process,  a  photon  transfers  all  of  its  energy  to  an 
electron,  which  is  "knocked  on"  and  dissipates  its  energy  in  the  immediate 
vicinity.  The  probability,  per  unit  length  of  travel,  that  a  photon  will 
experience  such  a  collision  is  denoted  (Ey  ).  When  the  atom  from 

which  the  electron  was  ejected  acquires  another  electron  in  its  place, 
an  x-ray  of  low  energy  is  emitted j  however,  about  95  per  cent  of  these 
are  reabsorbed  inside  the  atom,  and  an  "Auger"  electron  of  low  energy 
is  ejected  from  among  the  outer  valence  electrons.  This  electron  also 
dissipates  its  energy  in  the  vicinity.  Various  other  energy  "dribbling" 
processes  may  occur,  but  the  net  effect  is  that  all  of  the  energy  of  the 
primary  gamma  photon  is  dissipated  as  heat  in  the  vicnity  of  the  collision. 
Thus  the  fraction  of  initial  gamma  energy  Ey  dissipated  is  unity. 

Consequently,  the  heating  rate  due  to  photoelectric  collisions  of  gamma 
photons  is  given  by  ^  ^ 


where  E  ^  denotes  the  energy  of  the  gamma  photon,  and  ^y(E  y )  denotes 
the  energy  dependent  gamma  flvac. 
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Photons  having  energies  in  excess  of  1,02  mev  may  interact  vd.th  the 
coulomb  field  of  a  nuoleus  and  produce  an  "electron  pair". .  The  probabi¬ 
lity,  per  unit  length  of  travel,  of  pair  prodiiction  is  denoted  by  <^^pp(E  y) , 
The  resulting  beta  particles  dissipate  their  kinetic  energy  near  the 
point  of  production.  The  positron  ultimately  combines  with  an  electron 
forming  two  photons  having  energies  approximately  0,51  mev.  These  may 
penetrate  some  distance  from  the  point  where  the  pair  production  occurred. 
The  nucleus  Involved  in  the  collision  acquires  some  energy  also.  The  net 
result  is  that  (E^-  1.02)  mev  of  energy  is  dissipated  at  the  point  of 
collision.  Thus  ^^p>  tlie  fraction  of  photon  energy  Ey  dissipated  locally 
becomes 


^pp  -  (v  1-02) Ar  ;  1.02 


23 


The  heating  rate  due  to  pair  production  takes  the  form 

,00 


^5-pp  ‘^y/.'pp^y)  (E3')aK.y  24 


The  lower  limit  is  taken  as  zero  for  later  convenience,  (Ey)  is 

zero  for  Eyless  than  1.02, 

Compton  scattering  is  the  scattering  resulting  from  an  elastic  collision 


For  this  process,  Compton  derived  the  following  relation  between  the 
energies  E^  and  of  the  photon  before  and  after  collision  and  the 
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scattering  angle 


_ 1 _ 

ltE^(l-Cos^)/0.51 
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Since  the  collision  is  elastic,  kinetic  energy  is  conserved,  and 


E>ff  =  E^'-E^'  -  1  — S!: 


The  probability  that  a  photon  will  be  scattered  through  an  an^e  ^ 
into  angle  d-6-  is  given  by  the  KLein-Nlshina  formula  (  ^  ),  which 
may  be  put  into  the  form 

P(d)d&  -  KCE^  /^)2  -h  (^/V)  -  Si^SinPdO 


where  K  is  a  normalization  constant  such  that 


=  1 
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and  may  be  evaulated  by  substitutions  from  Equations  25  and  27  followed 
integration. 

The  fraction energy  E^of  the  incident  photon  that  appears 
as  kinetic  energy  E^  of  the  beta  particle  is  given  by 


^y-c  (%) 


(1  -  Ey/E2r)E^P(  9)d9 
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Substitution  from  Equations  25,  26  and  27  yields  an  integral  that  has 
not  been  put  into  a  convenient  algebraic  form  but  has  been  evaluated 
(See  Equation  32 et,  seq.). 

The  energy  deposition  due  to  Compton  scattering  may  be  formulated 


27 


12 


follows 


dE^r 


Summing  the  gaiimia  heating  contributions  given  by  Equations  29,  24  and  22, 


SVp^pe-fi^P/^Ppf^c^c)  Ey,j2fy(E^)  dE^ 


where  the  Ey  *s  have  been  omitted  in  the  parenthesis. 

Substituting  from  Ec. nation  23  and  rcarnfinglng  ^  but,  keeping  in  mind 
tha^,^pp  vanislies  for  energies  less  than  1.02  mcv.  and  that  ^ype 
is  unitj’-,  gives 


^/'pe^PP-f^A'c)  E^^^(%)dE^-  1.0^p<^,(Ey)dS 


The  quantity  (^«pe+y/^p  4- )^appearing  in  the  ],eft  hand  integral, 
divided  by  the  dcns.ity  of  the  medium,  has  been  graphed  hy  Snyder  and  Povrell 
(18)  as  a  function  of  Ey  for  veirious  materials.  Their  coefficient^- 
is  identical  with  the  quantity  in  the  parenthesis.  It  is  called  the  energy 
absorption  coefficient.  If  the  so-called  "annihilation  photons"^  resulting 
from  the  neutralization  cf  the  posltron^are  assumed  to  be  absorbed  in  the 
vicinity  of  the  point  of  pair  production,  one  may  neglect  the  right-hand 
integral  provided  the  annihilation  photons  are  not  included  in  the  gamma 
flux  )Ziy  (Ey  ). 

^  Denoted  elsewhere  by 


13 


MEUTRON  FLDX  DISTRIBUTION 

The  neutron  diffusion  and  transport  theories  yield  descriptions  of  the 
spectral  and  spatial  distributions  of  the  neutron  fluxes,  (^n)* 
theories  have  been  discussed  exhaustively  elsewhere  (  5  )  and  will  not  be 

considered  here  in  detail.  A  one  or  two-group  diffusion  calculation 
vdll  suffice  to  describe  the  "gross'*  variation  of  the  thermal  flux  in  many 
thermal  reactors  with  sufficient  precision  for  heat  generation  calculations. 
Let  ^(th,0)  be  the  flux  at  some  reference  point  in  a  homogeneous  reactor 
(say  the  center),  and  let  f  (  r  )  denote  a  function,  obtained  from  diffusion 
theory,  describing  the  spatial  variation  of  Ai  where  r  denotes  the 

radius  vector  from  the  reference  point,  such  that 


=  Pn  f(r) 

By  Eqi«itlon  3  -  it  follows  (using  the  notation  of  Equation  33)  that 


Gf  (e)  =  A  ^  (th,0)  f  (r)  (th,r) 


34 


where  A  is  a  proportionality  factor.  The  total  heat  output  from  the  reactor 
due  to  fission  fragments,  Qf,  is  obtained  by  integrating  Gf(r)  over  the  volume 
of  the  fuel-bearing  part  of  the  core. 


Of 


Gf  (r)  dV 


v/Vf 

where  Vf  denotes  the  volume  of  the  core. 

Substituting  from  Ecpiation  34  and  solving  for  A  j6n  (thjQ)  gives 


35 


A  (th,0) 


f  (2:)^f 


(th,r)  dV 


36 
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Putting  this  back  into  Equation  34  gives 


Gf  (e) 


f  (r)2.f 


yj 


The  energy  available  from  fission  ■^agments  (kinetic  and  beta  decay) 
is  about  175  mev*  per  fission.  The  total  energy  per  fission  is  about  200  mev. 
Thus  Qf  is  about  87  per  cent  of  the  total  reactor  heat  output.  The  bal¬ 
ance  of  is  generated  by  fission  gammas  (8  mev)^ fission  fragment  decay 
gammas  (6  mev),  fast  neutrons  (5  mev),  capture  gammas  and  decay  betas  (7  mev). 
The  particles  bearing  this  energy^about  25  mev)  tend  to  leak  out  of  the  fuel¬ 
bearing  core,  A  goodly  portion  of  their  energy  is  released  in  the  coolant, 
reflector,  thermal  shield,  pressure  shell,  and  biological  shield.  However, 
in  designing  the  fuel-bearing  core,  it  is  sufficient  and  conservative  to 
assume  that  all  of  the  reactor  heat,  ^  homogeneous  reactor  is  released 

in  the  fuel,  and  that  it  has  the  same  spatial  distribution  as  Qf. 

The  function  £  (r)  has  been  evaluated  by  Glasstone  and  Edlund  (  5  ) 
for  certain  bare,  homogeneous,  thermal  reactors  of  simple  shape,  with;^  f (th,V) 
uniform. 


Table  1 

Geometry 

f(r) 

Go/Gav 

Sphere  of  radius  R 

Sin  (7^'r/R)/(7rr/R) 

3.29 

Cube  of  side  A 

Cos  (7rx/4)  Cos  (7>y/A)  Cos(rrz/A) 

3.64 

Right  cylinder  of  radius  R 

Jo(2.405r/R)  Cos  (TT  z/2R) 

3.87 
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By  integrating  the  given  fiinctions  over  the  corresponding  geometric 
forms,  the  ratio,of  the  maximum  heat  release  rateii  G55.  to  thei.^terige  . 
rate,G^^jis  readily  chtained.  Values  are  listed  in  Table  1.  The  presence  of 
moderating  reflectors,  control  elements,  structural  materials,  coolants, 
ducts,  and  non-uniform  loading  (variable  )  modify  these  simple  distri¬ 
butions.  Space  does  not  permit  treatment  ofjthese  effects,  but  it  is  impor¬ 
tant  to  reduce  so  that  a  larger  fraction  of  the  core  may  operate 

at  or  near  the  limiting  conditions  and  bum-up  of  nuclear  fuel  and  accom¬ 
panying  effects  will  be  more  imiform, 

Glasstone  and  Edlund  point  out  that  the  functions  Jq(x)  and (Sin  x)/x 
are  rather  similar  to  Cos  x,  and  to  a  good  first  approximation  may  be 
replaced  by  Cos  x. 

In  order  to  increase  resonance  escape  in  low  enrichment  uranium  reactors, 
the  fuel  is  sometimes  segregated  from  the  moderator.  Fast  neutrons,  released 
by  fission  in  the  fuel,  have  a  high  probability  of  escaping  into  the  modera¬ 
tor  where  they  are  quickly  slowed  down  without  exposure  to  resonanct  absorp¬ 
tion  in  In  the  OHIIL  Graphite  Reactor,  this  result  is  achieved  by  using 

metallic  fuel  rods  about  1  inch, in  diameter  arranged  in  an  8  inch  square 
array  in  a  matrix  of  graphite.  The  fuel  has  negligible  moderating  proper¬ 
ties;  neutrons  are  "thermalized"  almost  entirely  in  the  graphite  and  dif¬ 
fuse  thence  back  into  the  fuel,  there  causing  further  fissions.  As  a  re¬ 
sult,  there  is  a  local  gradient  in  the  flux  density  both  in  the  moderator 
and  in  the  fuel  superimposed  on  the  gross  fluk  of  Eq.  33«  Some  reactors, 
which  are  heterogeneous  from  a  phase  standpoint,  such  as  the  MTR,  are 
really  homogeneous  from  a  nuclear  standpoint,  there  being  negligible  gradients 
in  the  neutron  flux  in  the  moderator  and  fuel. 


Neutron  diffusion  theory  may  be  employed  in  the  calculation  of  therm^ 
flux  distribution? in  fuel  elements  in  heterogeneous  reactors*  The  basic 
equation  is  (  5  ) 

(th,r)  -  ^  (th,r)  4-  S^(th,r)  =  3j/5t  3B 

where  n  is  the  number  of  thermal  neutrons  per  cco,  Sji  is  the  volume-source 
of  thermal  neutrons  (due  to  slowing  down  of  fast  neutrons)  in  neutrons  per 
cc.  per  second,  D  is  the  neutron  diffusion  coefficient,  Is  the  neutron 
absorption  coefficinet  (the  so-called  "macroscopic"  absorption  cross  section), 
and  t  denotes  time*  In  the  steady  state, "^n/^t  is  zero;  in  the  fuel,  is 
zero.  The  solution,  in  polar  coorda nates  (corresponding  to  the  case  of  an 
infinitely  long  fuel  rod  immersed  in  a  Tiniform  external  fl\ix),  is 

4  (th,r)  =  39 

where  t  is  the  distance  from  the  axis  of  the  fuel  rod  and  B  is  an  arbi¬ 
trary  constant  whose  magnitude  is  proportional  to  the  power  level  at  which 
the  reactor  is  operating  and  >diich  may  be  evaluated  by  multiplying  the  gross 
flux  distribution j  fC^)  of  Equation  35,  by  the  local  variation  {l^ji(th^)  and 
proceeding  through  ■ttie  operationa  indicated  in  Equations  34935  »  andh96  to 
evaluate  the  product  AB^tth^fi)  terms  of '  %  „ 

It  is  customary,  and  conservative  practice,  to  assume  that  nil  of  the 
fission  energy  (200  mev  per  fission)  is  released  in  the  fuel,  even  though 
the  fast  neutrons  and  a  significant  fraction  of  the  fission  gammas  and 
fission-fragment  decay  gammas  release  their  energy  elsewhere.  In  estimating 
heat  release  in  moderator,  coolants,  pressure  shells,  and  shield^  however. 

It  is  necessary  to  estimate  these  leakages  from  the  fuel  rather  precisely* 

A  fair  estimate  of  the  fast  neutron  spectral  distribution  may  be  had  from 


17 


the  Fermi  age  treatment  (  5)  of  simple  modelsi  e^go,  an  infinite,  hetero- 

235  238 

geneous  matrix  with  negligible  capture  of  fast  neutrons  in  TJ  and  U 
and  only  weak  captxxre  in  the  moderatoro  It  is  also  assumed  that  the  flux 
of  uncollided  neutrons  from  fission  is  uniform  throughout  the  moderatoro 

Alternatively,  the  Monte- Carlo  technique  (  17  )  permits  treatment 

of  the  spatial  variation  of  the  fast  flux  for  the  case  of  strong 
absorption  in  fuel,  moderator,  and  control  elements,  and  leakage  out  of 
the  reactor,  provided  the  fission  rate  density  in  the  fuel  can  be  esti¬ 
mated  by  some  means,  Sogo,  by  a  multi-group  diffusion  calculation  on  a 
"homogenized"  model  of  the  coreo  In  the  Monte-Carlo  method,  the  events 
in  a  neutron  "life  history"  are  selected  by  sampling  from  a  table  of  ran¬ 
dom  numbers*  Initial  neutron  energies  are  selected  from  the  fission  neutron 
spectral  distribution,  a  flight  direction  is  established,  the  coordinates  of 
the  first  collision  are  found,  kind  of  nucleus  and  type  of  collision  is 
selected,  post-^llision  neutron  energy,  and  new  flight  direction  are 
determined.  The  history  is  followed  until  the  neutron  is  absorbed  or 
leaks  out  of  the  reactor.  The  process  is  repeated  until  a  sufficiently 
large  number  of  cases  has  been  obtained  to  permit  a  statistical  ahaljralsj 
to  obtain  the  fast  flux  spatial  and  spectral  distribution. 

The  technique  is  ^plicable  to  the  calculation  of  fast  fl\uces  in 
moderator,  core  vessels,  shields,  reflectors,  thermal  shields,  pressure 
shells,  biological  shields,  etc.  Unless  simplifications  are  introduced, 
especially  geometric  symmetries,  an  excessive  number  of  histories  may  be 
required,  even  though  an  electronic  computer  is  used. 
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GAMMA  RAY  FLUX  DENSITY  DISTRIBUTIONS 

Compton  scattering  greatly  complicates  the  estimation  of  gamma 
flux  densities.  The  process  is  not  isotropic  (Equation  27),  and  the 
energy  decrement  of  the  scattered  photon  is  strongly  dependant  on  the 
scattering  angle  (Equation  2) .  Photons  scattered  through  small  angles 
retain  most  of  their  original  e  nergy  and  may  penetrate  t  o  remote  regions 
of  the  reactor.  Photons  scattered  through  larger  angles  suffer  large 
energy  losses,  and  their  penetrating  power  is  greatly  reduced. 

The  important  gamma  sources  are  the  primary  gamma  soiirces  (the 
fission  process,  the  decay  of  fission  fragments,  neutron  captures, 
inelastic  scattering  collisions  of  neutrons)  and  the  secondary  source 
(Compton  scattering  of  photons).  At  many  points  of  interest,  indeed, 
at  most,  the  contribution  to  the  gamma  flux  made  by"  scattered  photons 
is  many  times  greater  than  that  of  unscattered  photons. 

In  principle,  the  gamma  flux  may  be  constructed  by  integration 
over  the  gamma  sovirces.  Let  r  denote  the  vector  coordinates  of  a 
point  in  the  reactor,  and  let  f6^  (Ejr,  c)  denote  the  gamma  flux  den¬ 
sity  in  photon-cm»  per  second  per  cc,  per  mev.  Let  £*  denote  another 
point  in  the  reactor,  let  a  volume  increment  about  £*  be  denoted  by 
dV*,  and  let  the  gamma  source  strength  in  photons  per  second  per  cc. 
at  c’  be  denoted  by  Sy  (jf')®  Let  P(Ey  ,r')  be  the  conditional  pro¬ 
bability,  per  mev., that  a  photon  will  be  emitted  from^or  scattered 
out  of  dV  with  energy  E^-.  Let  P(£*,  z)  t>e  the  probability,  per 
steradian,  that  the  photon  will  be  emitted  in  the  direction  of  r. 

The  probability,  per  cm.  of  travel,  that  a  photon  will  be  removed 
from  the  beam,  denoted  simply  the  sum  of  the  proba^ 

bilities  of  the  photoelectric, pair  production,  and  Compton  scattering 


processes. 
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(E<r) 

The  probability,  P(R*-^  r),  that  a  photon  will  arrive  xmcoUided  at 
£  is  thus  simply  (  in  homogeneous,  isotropic  media) 

?(£'—>•£)  =  exp  ()  £*  “  X|)^  41 

where  /£*  -  £|  denotes  the  scalar  distance  between  £*  and  r.  It 
follows  that  ^ 

(Ey,  £)  =  / S^(r')P(Ey,E>)dVP(r»,E)  exp  [^(E^)  (l£’-£^42 

The  primary  gamma  sources  are  isotropic,  and  for  these  P(e'>£) 
has  the  value  l/^TTper  steradian.  In  the  case  of  the  secondary  gamma 
source,  P(3?',e)  is  a  complicated  function  of  the  "directed"  gamma  flux 
(Ey  r’),  which  is  defined  as  the  gamma  flux  at  e'  contributed 
by  photons  of  energy  E  moving  in  the  direction-^  «  It  has  tmits  of 
neutron-cm.  per  second  per  cc.  per  mev  per  steradian.  Thus  the 
construction  of  the  flux  (Ey  ,  £)  at  r  requires  a  knowledge  of  the 
more  complicated  flux  ,  e')  everywhere  else.  Consequently, 

a  general  formulation  of  e)  does  not  seem  feasible. 

The  Monte- Carlo  technique  may  be  applied  to  this  problem  (  14  )• 
Suppose  the  histories  of  sufficient  primary  photons  of  energy  E'  a^e 
determined  to  define  F  (E‘5  the  average  fraction  of  the 

energy  E«  of  photons  originating  af  e'  that  is  dissipated  at  e 
per  cc.  The  heating  rate  is  now  readily  formulated  as  follows. 


Spy(r’)dV'  P(EV)  dEf  E^  F(E^je'-^£) 


(v  /e/-  0 


i 


where  Sp  ^  (r * )  Is  the  primary  gamma  source  strength. 

In  principle,  r)  may  be  evaluated  for  any  case  howsoever 

complicated  in  geometry  or  heterogeneity,  but  the  sheer  magnitude  of  the 
computational  labor  nay  render  the  method  impractical.  The  simplest 
case  is  that  of  an  infinite,  homogeneous,  isotropic  medium.  Life 
histories  are  generated  in  the  usual  manner  and  an  account  is  kept  of 
the  radial  distance  of  each  collision  from  the  point  of  origin  and 
the  amount  of  energy  deposited.  From  these  data  it  is  easy  to  calcu¬ 
late  the  fraction  f(E',R)  dR  of  the  energy  E'  of  the  original  photons 
that  is  released  inside  a  spherical  annulus  of  radius  R  and  thickness 
dR.  The  vol\nne  of  such  an  annulus  is  47^^^  dR,  from  which  it  follows 
that  F(E*jr'-^r)  is  given  by 

F(E'iE--.E)  =  -  Elfl  u. 

47r  (  In’  - 

A  brilliant,  analytical  solution  of  the  foregoing  case  has  been 
obtained  by  Goldstein  and  Wilkins  (  6  )  along  lines  laid  down  by- 
Spenser  wriH  Fano  (  19  )•  Briefly,  the  method  consists  in  expanding 
the  directed  flux,  (Ey  ,  it  ,£),  in  a  series  the  terms  of  which  are 
products  of  Laguerre  and  Legendre  polynomials.  The  directed  flux  is 
integrated  with  respect  to  all  directions  to  remo-ve  the  angular  depen¬ 
dency,  multiplied  by  the  distance  from  the  origin  raised  to  integral 
powers,  and  the  product  is  integrated  over  all  space  to  eliminate  the 
spatial  dependency.  The  integration  yields  the  so-called  ”moments”  of 
the  flux,  which  are  used  to  evaluate  the  arbitrary  coefficients  in  the 
original  series  in  the  usual  way  by  application  of  the  boundary 
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conditions  and  utilization  of  the  orthogonal  properties  of  the  Legendre 
polynomial.  Having  constructed  the  flux  £)  due  to  a  monoenergetic, 

unit,  point  source  (one  photon  per  second  of  energy  E’)  at  £*,  Goldstein 
then  computes  the  fraction  of  the  energy  E'  released  at  Eo  This,  of 
course,  is  identical  with  F(E'^,  E»  — ^  r)  defined  above.  The  resiilts 
are  reported  in  terms  of  an  "energy  absorption  build  up  factor^ 

®a  where  E^  is  E'  in  the  notation  used  here^y^o  is^^(E'), 

r  is  ^|r'  —  E/J,  and  where  is  defined  in  such  a  way  that 

F(E',£'^r)  =  45 

where^  is  the  energy  deposition  coefficient  defined  in  Eqxiation  32, 

Values  of  ^  function  of  and^E'  -  £|  are  tabulated  for  a  variety 

of  materials  by  Goldstein  (  6  )  „  Taylor  (  20  )  has  proposed  an  empiri¬ 
cal  correlation  exponential  in  r*  -  e  which  is  conveniently  applied  to 
problems  in  gamma  heating  in  pressure  shells  eind  thermal  shields  (  1  ) , 

Conceivably,  the  function  F(E',|^'  -  r^oould  be  evaluated  experi¬ 
mentally,  Small  radiation  sources  of  various  energies  could  be  moved 
about  in  a  mock-up  of  the  system  under  study  and  the  heat  liberation 
at  various  points  of  interest  measured.  If  a  direct  thermal  measurement 
is  made,  the  measuring  instrument  would  have  to  be  extremely  sensitive 
or  else  the  radiation  source  strength  would  have  to  be  inconveniently 
large.  If  an  electronic  measurement  is  made,  the  spectral  distribution 
of  the  gamma  flux  would  have  to  be  measured  in  order  that  the  heating 
might  be  computed  by  Equation  43, 

The  mathematical  difficulties  in  estimating  gamma  heating  may  be 
considerably  reduced,  albeit  at  an  unknown  cost  in  accuracy,  if  one  is 
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■willing  to  utilize  the  straight-ahead  scattering  approximation,  investigated 
carefully  by  Huruitz,  et  al.  (ll).  Hie  assumption  is  made  that  photons 
suffer  energy  degradations  according  to  Equations  25  and  27,  but  that  the 
scattered  photons  fall  into  one  of  two  classes:  (a)  photons  scattered 
through  such  small  angles  that  the  de-vlations  from  the  line  of  flight  may 
be  ignored,  and  (b)  photons  scattered  through  such  large  angles  that  the 
resulting  low-energy  photons  have  little  penetrating  power  (high  pt)  and 
are  absorbed  near  the  point  where  they  suffered  the  scattering.  These 
assumptions  depend  upon  -the  following  arguments:  (a)  the  energy  distribution 
of  the  scattered  photons  shows  a  strong  forward  component,  (b)  the  angular 
distribution  also  shows  a  strong  forward  component,  (c)  distances  of  Interest 
in  heat  generation  calculations  are  measured  in  relatively  few  attenuation 
lengths  from  the  primary  source,  and  the  probability  that  a  photon  will 
suffer  many  scatterings  is  small,  (d)  photoelectric  absorption  Increases 
strongly  with  decreasing  photon  energy.  For  example,  of  the  photons 
scattered  out  a  beam  of  6  mev  photons,  those  scattered  through  angles  less 
than  50  degress  carry  90  per  cent  of  the  energy  of  all  the  scattered  photons. 

The  approximation  is  applied  by  breaking  the  gamma  soinrce  spectrum 
into  a  number  of  energy  groups.  The  group  having  the  highest  energy  Ej  is 
attenuated  betv/een  r’  and  r  by  the  simple  exponental  factor 
exp  -  £l)  ]•  The  next  group  loses  photons  in  proportion 

to  gains  photons  from  group  1  by  Compton  scattering.  Group  3 

receives  photons  from  both  groups  1  and  2,  and  so  on.  The  transfer 
of  photons  between  groups  is  a  function  of  the  energy-widths  of  the 
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groups  and  their  separation  in  energy,  and  is  computed  from  the  Klein- 
Nishina  relation  (Equation  27).  The  differential  equations  expressing 
the  photon  •’economy"  for  each  group  are  integrated  successively,  giving 
the  flux  at  r. 

Hurwitz,  using  this  multi-group  technique,  has  calculated  the  energy 

release  ar  r  diie  to  an  energy-distributed  source  at  r'  having  a  typical 

spectrum  and  found  that  the  results  could  be  expressed  in  terms  of  an 

average  "effective"  gamma  absorption  coefficient, (E^),  obtained 

by  averaging/<(!^ff  over  the  spectrum  of  the  source 

/ 

where  P(E')  is  the  spectral  distribution  and  where 

47 

is  a  term  added  to  the  energy  deposition  coefficien^/^  to 
account  for  the  energy  release  due  to  gamma  photons  scattered  to 
energies  below  0.25  mev,  which  are  assumed  to  be  absorbed  in  the  near 
vicinity  of  their  scattering.  It  is  evaluated  ty  integrating  the 
scattered  energy  ratio;  E^/Eygiven  by  Equation  25  with  respect  to  the 
angular  probability  given  by  Equation  27  as  follows: 

Ae(E«=  0.25) 

(E)  =  I  (^/E)P(  -Q)  48 


arrwitz's  results  may  be  put  in  the  form,  using  the  present  notation, 

ojt)  =  /  Sp,r(S')  dV  49 

jv  4Tr(ir'  -  eF  t 
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PRIMARY  GAMMA  SOIffiCES 


The  fission  gamma  and  fission  fragment  decay  gamma  sources  are  both 

proportional  to  the  fission  density  in  stationary  fuel  reactors, 
pOO 

^  dE^  50 

J  0 

where  N^fis  the  number  of  photons,  fission  or  decay  as  the  case  may 
\5e  emitted  per  fission.  Gamble  (  4  )  reports  that  7,5  prompt  gamma 

photons  are  released  per  fission.,  that  these  have  an  average  energy  of 
about  1  mev,  each,  and  that  their  spectral  distribution  is  approximately 
exponential, 

P(E^)  =  exp  (  -E^)  51 

This  relation  does  not  hold  at  the  extremes  of  the  energy  range;  however, 
the  contribution  from  photons  with  energy  above  ten  mev*  is  negligible, 
while  gammas  emitted  with  energies  less  than  0,1  mev<  carry  less  than 
one  per  cent  of  the  total  energy o 

The  fission  fragment  decay  gamma  spectrum  may  be  constructed  from 
a  knowledge  of  the  fission  yields  of  the  various  nuclides  and  their  decay 
schemes,  Moteff  (  13  )  has  tabulated  these  data  for  some  nuclides  which 
have  appreciable  fission  yields,  which  emit  gammas  with  energies  in  excess 
of  0.04  mev,,  and  which  have  half-lifes  in  excess  of  30  seconds.  Several 
other  radioactive  fission  products  are  known,  and  it  seems  likely  that 
considerable  gamma  energy  is  radiated  by  nuclides  having  half-lifes 
lying  between  Gamble’s  "prompt"  period  (^Ji.  microseconds)  and  those 
which  have  been  observed.  In  all,  about  6  mev.  is  released  per  f±saiion(2l). 

The  situation  is  fairly  satisfactory  with  respect  to  data  on  the 
spectrxim  of  the  capture  gamma  rays  emitted  by  a  nucleus  in  returning 


to  the  ground  state  following  absorption  of  a  neutron.  Mittleman  (  12  ) 
has  tabulated  the  information  available  as  of  October,  1953.  The  strength 
of  this  gamma  source  is  proportional  to  the  local  neutron  flux  and  to  the 
capture  coefficient  ^c.  Equation  50  may  be  used  by  replacing^^  ^2%“ 

Hardly  any  data  are  available  on  the  spectrum  of  gamma  rays  emitted  fol¬ 
lowing  the  fast  neutron  inelastic  collisions.  Grace,  et  §1 . (  V  ) 
report  some  measurements  on  the  scattering  of  2o5  mev  neutron^p 
HLizard  (  3  )  states  that  it  may  be  possible  to  predict  these  spectra 
in  the  same  way  that  spectra  of  ordinary  gamma  emitters  are  predicted 
(  8  ),  but  nothing  seems  to  have  been  published  on  this. 

CONCLUSION 

The  local  rates  of  heat  relea.se  in  nuclear  reactors  have  been  for¬ 
mulated  in  terms  of  the  local  neutron  and  gamma  fluxes.  Certain  special 
methods  of  estimating  these  fluxes  were  indicated  and  their  application 
to  the  heating  problem  discussed  briefly. 
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